Electro-viscous effects in capillary filling of nanochannels 
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Abstract 

We theoretically examine the widespread hypothesis of an electro-viscous origin of the increase 
in apparent viscosity observed in recent experiments on capillary filling of nanochannels. Including 
Debye-layer corrections to the hydraulic resistance we find that the apparent viscosity reaches 
a maximum in the mesoscopic regime where the channel height (or more generally the hydraulic 
radius) is comparable to the screening length. However, for realistic estimates of central parameters, 
we find that the electro- viscous contribution to the apparent viscosity is at most a 1% effect. 
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Passive nanofluidic devices have recently been explored in a number of contexts aiming 
at both fundamental studies such as DNA stretching; as well as for applica- 

tions in optofluidic dye lasers {(J. During the past couple of years this has also stimulated 
basic experimental investigations of fluid dynamics particularly addressing possible new phe- 
nomena associated with the cross-over from a macroscopic (i.e. micrometer scale) to the 
nano regime. The list relevant length scales and phenomena associated with this cross over 
includes the Debye screening length A^ (typically of the order 10 nm), the zeta surface 
potential ( (typically sub-fe^T and of the order meV), and the issue of slip versus no-slip on 
the inter- molecular sub-nanometer length scale 8 Q. 

Of particular importance is to achieve a more complete understanding of fluid dynamics 
in the mesoscopic regime, where relevant length scales of the fluid, e.g. the Debye screening 
length, is comparable to device dimensions, for example the channel height h, while still 
exceeding the molecular length scale, i.e. 5 <C \d < h. This mesoscopic regime is typically 
entered when capillaries with cross sectional dimensions of approximate 100 nm and below, 
is infilled with buffer solutions with an ionic strength relevant for carrying typical biological 
samples, e.g. DNA. A typical salinity in the mM range results in Ad ~ 10 nm Q. 

This has motivated several basic nano-scale experimental investigations on the classical 
problem of capillary filling, see Fig. [TJ thus testing the classical Washburn formula [8] 

9 . . 7/1 cos # 

At) = (i) 

which is known to provide a well-established description of the filling dynamics in capillaries 
with cross-sectional dimensions of mm size and down below 1 /mi. Here, x is the length 
of advancement along the channel of the liquid front, 7 is the water- air surface tension, 6 
is the contact angle, and r] is the viscosity. The Washburn formula is derived within the 
framework of the Stokes equation for a charge-neutral liquid where no surface-charging and 
electro-osmotic phenomena occur. Eq. (pQ) implicitly relies on the continuum hypothesis 
which limits its validity for fluid dynamics at the inter-molecular length scale while it is still 
expected to work in the mesoscopic regime where physical mean-fields remain well-defined 
and continuum formalism applies. 

Experimental observations support an unexpected increase in the filling time, beyond the 
predictions of Eq. (jTJ, in shallow nanoslit devices for slit heights of approximate 100 nm 
and below 0, [lOj. In all reported experiments, a systematic deviation from the behavior 
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predicted by the Washburn formula is observed: Although the observed dynamics maintain 

the proportionality between x 2 and t, the constant of proportionality a does not scale linearly 

with the channel height h for slit channel heights of approximate 100 nm and below. 

In absence of a detailed theoretical model, different possible interpretations have been 

presented for the unexpected height dependence of a(h). Tas et al. speculate that the 

phenomena is of an electro-viscous origin 91 while Thamdrup et al. show that in their 

□ 

experiments the phenomena correlates with the formation of air bubbles [10jj. Both inter- 
pretations implicitly assume that changes in hydraulic resistance -Rhyd play a central role. 
For this reason, the changes in observed filling time are naturally presented in terms of an 
increase in apparent viscosity rj, entering the Washburn formula. 

In this Letter we work out the electro-viscous correction to the apparent viscosity r? a , p P 



and discuss the result and consequences in the context of the experiments in Refs. [9 



io|. 



Our theoretical account of the e 
to the hydraulic resistance [ll|, 



ectro- viscous correction includes the Debye-layer correction 
1^ | . When decreasing the height h the apparent viscosity 



increases in qualitative agreement with the experiments, reaching a maximum when the 
height of the channel h is comparable to the Debye screening length Ad, and then drops 
off again when further reducing the height beyond the Debye screening length, see Fi g, [p 
This picture agrees with earlier theoretical work on retarding electrokinetic counter flow [13 ]. 
However, when scaling the experimental results with the Debye screening length the quanti- 
tative agreement with experiments is poor and the prediction of a maximum is not consistent 

nn 

with experiments [9|, IIOJ • This suggest that there are other phenomena contributing to the 
apparent viscosity, thus supporting the explanation by the formation of air bnbb.es Q over 
the electro- viscous hypothesis ]9J. 

The basic physics behind the electro- viscous hypothesis has already been discussed by Tas 
et al. jo |, but here we repeat the arguments for completeness. When water and other elec- 
trolytes flow in a channel, chemically induced charge-transfer takes place at the channel walls 
rendering a charged Debye layer of width Ad in the liquid near the wall, and a compensating 
oppositely charged wall. This electric phenomena combines with the purely viscous flow and 
results in electro- viscous effects. In the spirit of the derivation of the Washburn formula , 
the apparent viscosity of course stems from an apparent increase in hydraulic resistance i?h y d 
due to an electro-osmotic counter flow. Since the hydraulic resistance scales linearly with 
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the viscosity we have the following for the Debye-layer correction to the apparent viscosity 

Vzpp = -Rhyd(C) (2) 

v RhydiC = o) 1 ' 

where ( is the zeta-potential at the surface of the channel. 

The corrections from a finite zeta potential follows from our recent general linear-response 
theory of combined mass and charge transport in long straight channels III, |l2j with the 
electrolyte being subject to the application of an external pressure drop Ap and voltage 
drop AV along the channel. These potential drops results in a liquid flow rate Q and an 
electric current / given by the 2x2 conductance matrix G as 



(3) 



with the matrix being subject to an Onsager relation, G 12 = G^i- in the absence of a zeta 
potential we of course have Gi 2 (( = 0) = G 2 i(C = 0) = so that Eq. (T5J) becomes 

-i -l 
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(4) 
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where we have used that i?h y d(C) — -Rn(C) with the resistance matrix R = G^ 1 given by 
the inverted conductance matrix. In applying the result in Eq. ([3]) to the present problem 
of capillary filling we have implicitly utilized that the time scale for charge-transfer and 
charging of the Debye layer is several orders of magnitude smaller than the time-scale for 
moving the meniscus forward the distance of one Debye length. This difference in time 
scales secures that at any instant the charge distribution has reached its equilibrium, and 
consequently the electric current is zero, 1 = 0. 

We now proceed by calculating the conductance matrix elements by analytical means. 
For h <C w the problem becomes quasi one-dimensional and within the Debye-Hiickel limit 
it can be solved exactly. In particular we have 

h 3 w 
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where L is the length of the channel, a = eD/\ 2 D is the conductivity of the liquid, D is the 
diffusion constant of the ions, and e the dielectric function of the liquid. Long, but tedious 
manipulations then leads to 
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v = *> - ^ - 'If*** (7\ 
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which is correct to lowest order in the the dimensionless parameter 

6 2 C 2 _ < 2 _ wl 

r]a\ 2 D r]D eD 

expressing the square of the ratio between the electro-osmotic mobility and the viscous 
mobility The latter form of x is perhaps most appropriate from an experimental point of 
view, as the electro-osmotic mobility fi eo = eC/v can b e determined by measurements easier 
than the zeta potential ( itself. 

In Fig. [2] we show the apparent viscosity for increasing values of \. Most importantly it 
should be emphasized that the electro- viscous correction is maximal for h ~ 3.212Ad which 
is a value that can be found numerically. This is the mesoscopic regime where the Debye 
screening length A^ becomes comparable to the channel height h, while still exceeding the 
inter-molecular length scale 5. 

The p.„t also iae.uoes exposal data reproduced f r0m Ref, 0, Q for both a CM 

n n n 

NaCl solution with ~ lnm |9|, HJ] as well as for demi water with A^ ~ 200 nm [9(. 
While the experiments with a NaCl solution belongs to the sup-mesoscopic regime with 
h ^> Ad the demi- water experiments are in the opposite sub-mesoscopic limit with h < Xd- 
As seen, the demi-water results are not consistent with the hypothesis of a correction of an 
electro-vicious origin as derived above. We note that this conclusion is quite solid, since the 
demi-water data would belong to the sub-mesoscopic regime (to the left of the peak) even 
if we for some reason would allow for an one-order of magnitude reduction in the estimate 
of the screening length, i.e. A^ ~ 200 nm — > 20 nm. The dashed line shows a fit of Eq. (J6]) 
to the NaCl data with % ~ 16 used as fitting parameter. However, we emphasize that this 
high value of \ is completely unrealistic and in fact for a zeta potential in the meV range we 
find that x ^ 10~ 3 thus leaving the estimated relative electro-viscous correction a sub-1% 
effect at most. 



Nanofluidic capillaries have recently Deen usea m optonuiaics [yj ana m applications lor 
stretching of DNA molecules The majority of these experiments employ 



3een used in optofluidics [6] and in applications for 



nanochannels with roughly a square cross section. The above quasi ID results may be 
generalized to more arbitrary 2D channel cross sections with an aspect ratio of order unity. 
Following Ref. [l2| we straightforwardly arrive at the approximate result 



2 / R 



, 8 ^ V> 7 r x (8) 

where is the modified Bessel function of the first kind of order n and R = 2A/P is 
the hydraulic radius with A being the cross-sectional area and P being the length of the 
perimeter. For a circular cross section with radius a we have R = a and Eq. (jSJ) is exact in the 
limit of a small x- Likewise, for a square channel with side length a we get R = a/2. Eq. (jSJ) 
gives rise to the same qualitative dependence on R/Xd as in Fig. [2] and our above estimates 
and conclusions remain unchanged for nano channels with a true 2D electro-hydrodynamic 
flow profile. 

While we are confident that the electro- viscous effect would reveal itself in a clean exper- 
iment (i.e. in the absence of other competing effects) with a true nano-scale slit (h ~ 3Ac) 
the strength of the effect is far too weak to account for the experimental observations. This 
strongly supports the existence of other dominating effects, as also observed in recent studies 
where the apparent viscosity was found to correlate with the presence of air bubbles which 
tend to increase the hydraulic resistance [lo| . 

In conclusion, our calculation of the Debye-layer correction shows that for realistic es- 
timates of central parameters, the electro-viscous contribution to the apparent viscosity is 
below a 1% level. 
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FIG. 1: (Color online) Schematic illustration of a nanoslit-capillary filling experiment. A nanoslit 
capillary, of height h < 300 nm and width w several (im is filled with liquid from an attached 
reservoir by capillary action. The length of the fluid plug - the position of the liquid meniscus - 
x(t) is recorded as function of time t. 
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FIG. 2: (Color online) Plot of the apparent viscosity r] app /r] versus normalized channel height h/Xu 
for increasing values of the the ratio between the electro-osmotic mobility and the viscous mobility 
X = 0.1,0.5,1, and 2 (solid lines from below). The solid data points show experimental results 
reproduced from Ref. [9] for demi water with Xd ~ 200 nm (left group of data) and 0.1M NaCl 
solution with Xd — lnm (right group of data). The open data points show corresponding results 
for 0.1M NaCl reproduced from Ref. 14] with the dashed line being a fit to the data with x ~ 16 
used as fitting parameter. 
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